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ABSTRACT 



Aims. Interplanetary Lyman a line profiles are derived from the SWAN H cell data measurements. The measurements cover a 6-year period 
from solar minimum (1996) to after the solar maximum of 2001. This allows us to study the variations of the line profiles with solar activity. 
Methods. These line profiles were used to derive line shifts and line widths in the interplanetary medium for various angles of the LOS with 
the interstellar flow direction. The SWAN data results were then compared to an interplanetary background upwind spectrum obtained by 
O STIS/HST in March 200 1 . 

Results. We find that the LOS upwind velocity associated with the mean line shift of the IP Lyman a line varies from 25.7 km/s to 21.4 
km/s from solar minimum to solar maximum. Most of this change is linked with variations in the radiation pressure. LOS kinetic temperatures 
derived from IP line widths do not vary monotonically with the upwind angle of the LOS. This is not compatible with calculations of IP line 
O profiles based on hot model distributions of interplanetary hydrogen. We also find that the line profiles get narrower during solar maximum. 

Conclusions. The results obtained on the line widths (LOS temperature) show that the IP line is composed of two components scattered by two 
^ hydrogen populations with different bulk velocities and temperature. This is a clear signature of the heliospheric interface on the line profiles 
. . seen at 1 AU from the sun. 
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1. Introduction et al (1995, 1998) report measurements made by the 

GHRS (Goddard High-Resolution Spectrograph) instrument 
The interplanetary background has been used to study the in- j^^j^ble Space Telescope. More recently, new spectra were 
terplanetary hydrogen distribution since its discovery in the obtained from the STIS intrument that replaced GHRS on HST. 
late 1960's (Thomas and Krassa 1970; Bertaux and Blamont 

1970). Lists of previous experimental studies of the interplan- S^'^'lying the IP Lyman a line profiles has proved to be of 
etary Lyman a background can be found in Ajello et al. (1987, ^"^^1"^^^ ^^'^'lying the heliospheric interface because eff^ects 
1994) or Quemerais et al. (1994). imprinted on the hydrogen distribution at a large distance from 



the sun are still visible on Lyman a line profiles seen at 1 AU. 
This is due to the fact that the flow of hydrogen in the helio- 
sphere is collisionless. 



The SWAN instrument on the SOHO spacecraft (Bertaux et 
al., 1995) is able to study the IP Lyman a line profile using hy- 
drogen absorption cells. These devices are used to scan the line 
profile taking advantage of the variation of the Doppler shift Izmodenov et al. (2001) have presented model computa- 
between flie cells and flie interplanetary hydrogen as the space- tions that predict how the hydrogen distribution is aff^ected by 
craft rotates on its orbit around the sun. Previous studies of the the heliospheric interface. The heliospheric interface is the 
IP line profile with H cells were made by the Mars 6 space- boundary region where the expanding solar plasma interacts 
craft (Bertaux et al., 1976) the Prognoz 5/6 UV photometers with the interstellar plasma. Using the derived hydrogen distri- 
(Bertaux et al. 1985) and the ALAE instrument on ATLAS-1 butions, Quemerais and Izmodenov (2002) calculated line pro- 
shuttle mission (Quemerais et al. 1994). hl^s including a complete model of radiative transfer effects. 

IP Lyman a line profiles can also be obtained from The derived line shifts and line widths were clearly modified 

high-resolution UV spectrometers in die Earth's orbit. Clarke heliospheric interface. 

In this paper, we continue the work presented by Quemerais 

Send ofprint requests to: E. Quemerais et al. (1999) about the line profile reconstruction technique 
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based on the SWAN H cell data. This early work was concen- 
trated on the first year of data. Here we show results obtained 
from 6 one-year orbits of SOHO around the sun. Some of the 
details of the data processing are not repeated here and the in- 
terested reader should refer to the original paper. Previous anal- 
yses of the SWAN H cell data were pubhshed by Costa et al. 
(1999), Lallement et al. (2005), and Koutroumpa et al. (2005). 

The first section briefly presents the line profile reconstruc- 
tion technique and the SWAN data. The following sections give 
the results obtained for the Une shifts and the line widths, as 
well as the variations during the solar cycle. The last section 
presents a STIS/HST upwind fine profile obtained in March 
2001 and compares the derived line shift and Une width with 
the SWAN results. 

2. SWAN H Cell Data from 1996 to 2003 

The SOHO spacecraft was launched in December 1995 from 
Cape Canaveral. The SWAN instrument started operation in 
January 1996 and has been active since then except for a few 
periods of time (June to October 1998 and January 1999). 

Characteristics of the SWAN instrument are given in 
Bertaux et al. (1995). The SWAN instrument is a UV photome- 
ter with a passband between 110 nm and 160 nm. The instru- 
ment is made of two identical units placed on opposite sides of 
the SOHO spacecraft (+Z and -Z sides). Each unit is equipped 
with a periscopic scanning mechanism that allows to point the 
field of view in any direction of the half sky facing the side 
the unit is attached on. The instantaneous field of view is 5° by 
5° divided in 25 pixels. Each pixel has a 1° by 1° field of view. 
Measurements are performed every 15 seconds, with at least 13 
seconds of integration time to get a good signal-to-noise ratio. 

Each unit is equipped with a hydrogen absorption cell. 
This cell is placed on the photon path to the detector. The cell 
is filled with molecular hydrogen and has MgF2 windows. A 
tungsten filament passes through the cell. When a current goes 
through the filament, H2 is partially dissociated into atomic hy- 
drogen, which creates a small cloud that can absorb Lyman 
a photons. Typical values for the optical thickness of the ac- 
tive cell are around 3 to 5. Descriptions of the observing pro- 
grammes and of the various subsets of data are detailed by 
Bertaux et al. (1997). Here we concentrate on the hydrogen 
absorption cell measurements obtained between June 1996 and 
June 2003. 

The most common observation programme of SWAN is the 
so-called full-sky observation, during which each sensor unit 
covers the complete hemisphere that is on its side by moving a 
two-mirror periscope mechanism. One full-sky observation is 
performed in one day. The data obtained by both sensors are 
then combined into one image of the whole sky at Lyman a. It 
must be noted that the areas of the sky viewed by each sensor 
overlap, which enables us to compare both sensors on a regular 
basis. SWAN performs these observations 4 times per week. 
During the first year, one of the four full-sky observations made 
each week was made with cyclic activation of the hydrogen 
cells. For these observations, we then obtained a full-sky image 
at Lyman a as well as a full-sky image of the reduction factor, 
which is defined below. The mechanism is kept fixed during 



both measurements, cell OFF and cell ON, before moving to 
another direction in the sky. 

The reduction factor R used in what follows, is a dimen- 
sionless quantity. It is the ratio of intensities measured in a 
given direction when the cell is on (absorption) and when it 
is off (no absorption). If we consider an incoming intensity 
expressed as I(Ac) where Ac is the wavelength in the cell rest 
frame and if we note T(Ac) the transmission function inside the 
cell, the reduction factor is defined by 



Ion /r -^Wc) T(Ac) dAc 



/_ J I{Ac) dAc 



(1) 



Note that the quantity A = I - R measures the absorbed 
fraction of the incoming intensity. The transmission function 
inside the cell can be approximated with excellent accuracy for 
low optical thicknesses (r < 10) by 

7'(aj = exp(-Te-^), (2) 
where the variable x is defined as the normaUzed frequency by 



Vc - 



X = 
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Vo l Ac-Ao \ 
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(3) 



The V variable represents frequency, A wavelength. Here, Ao is 
the wavelength of the Lyman a transition at 1215.66 A. The 
Doppler width of the cell, A/t^, is related to the temperature 
and the thermal velocity in the cell by 



A 1 , Ao 

l^Ad Vthcell = — 

c c \ Mh 



2kT, 



cell 



(4) 



The absorbing power of a hydrogen cell can be characterised 
by a quantity called the equivalent width, Wa, in wavelength 
units. 



f 



Ll - T(Ac)] dAc 



(5) 



Typical maps of reduction factor values are shown in Bertaux 
etal. (1997). 

Quemerais et al. (1999) have shown how H cell measure- 
ments over a 1-year period can be used to reconstruct the line 
profiles of the IP glow. The same technique was applied here, 
and we refer readers to Sect. 3 of Quemerais et al. (1999). We 
use the same terms and notations as in this paper. 

The data used in this work cover a much longer period of 
time than our previous work which was concentrated on the 
first year of SWAN data. To reconstruct the line profiles, we 
need a full rotation of the Earth on its orbit. For each year we 
selected data from early June of the year until end of May of 
the following year For instance, the spectra labelled 2000 were 
derived from measurements starting in June 2000 and ending 
in May 2001 . The reason for this sampling was the lack of data 
between June 1998 and October 1998. As a consequence the 
year 1998 (June 1998 to May 1999) is missing. 

As done by Quemerais et al. (1999), we use Eq. (5) to caU- 
brate the absorbing power of the cell. Hydrogen cells are known 
to age with time. This may be due to ageing of the tungsten fil- 
ament itself or due to trapping of hydrogen by the glass or even 
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to some leak. This results in a decreasing optical thickness of 
the cell for a given filament current level. Over a 6-year period, 
both H cells have aged in very different ways. The H cell in the 
unit attached to the -Z side of the spacecraft seems to have lost 
all absorbing power very rapidly in 1999. This suggests that 
the cell has lost its H2 rather rapidly. The H cell in the other 
unit (+Z side) still retains most of its H2 cloud as shown by 
the strong absorption still seen in the 2005 data. However, the 
optical thickness has decreased with time. This was calibrated 
by determining the equivalent width of the cell (Eq. 5) from the 
data (see Sect. 4.1 of Quemerais et al. 1999). 

The values for the equivalent width of the +Z hydrogen ab- 
sorption cell are given in Tabled Note that we use the transfor- 
mation into velocity units given by Eq. (1 1) of Quemerais et al. 
(1999). These values were used to determine the line profiles 
as presented in section 4.2 of Quemerais et al. (1999). 

3. Line-of-sight velocities 

This section presents the results found for each of the 6 orbits 
analysed. As mentioned earlier, each data set starts in June and 
ends at the end of May of the following year. There is roughly 
one map per week, which means that about 50 maps are used 
to derive a line profile in each direction of the sky. In 1996 
and 1997, both H cells were active, so we have a full-sky ve- 
locity map for both orbits. However, due to the leak in the -Z 
unit H cell, only the northern ecliptic hemisphere is available 
from 1999 to 2003. We do not show individual plots of the line 
profiles. Examples are given in Fig. 5 to 7 of Quemerais et al. 
(1999). 

The line-of-sight (LOS hereafter) velocities shown here 
correspond to the mean Doppler shifts of the line profile ex- 
pressed in terms of velocity in the solar rest frame. If the line is 
not symmetrical, as in Quemerais and Izmodenov (2002) for 
example, then there is a small difference between the mean 
Doppler shift and the maximum of the line. 

Following Quemerais et al. (1999), we express the relation 
between velocity v projected on the LOS and the wavelength in 
the solar rest frame as 

i-^-^(y-[/) = ^(^) (6) 

where v is the projection of the atom velocity V on the line 
of sight U . Note that the direction of the LOS is opposite the 
direction of propagation of the photon. In that case, the LOS 
velocity < v > is given by 

1 c 

< v>= V I(v) dv ^ — < A- Ao> . (7) 

toff J-00 -io 

3.1. Comparison of the first two orbits 

Quemerais et al. (1999) show the velocity map derived from 
data obtained in 1996 and early 1997. This map corresponds to 
the minimum of activity of cycle 22. In the case of a flow with 
constant velocity, the velocity projected on the LOS is simply 
Voo ■ cosO, where 6 is the angle of the LOS with the upwind 




Fig. 1. full-sky map of the velocity difference between the 
1996 map and the 1997 map. The data are shown in ecliptic 
coordinates (west ecliptic longitude). The blue color in the up- 
wind side shows that the upwind velocity in 1997 is slightly 
slower (in modulus), by 0.4 km/s. 

direction and is the constant bulk velocity. The actual ve- 
locity maps observed by SWAN are more complex than this 
simple case. 

FigureQldisplays the difference between the velocity maps 
of 1996 and 1997. The variation in the upwind velocity is a de- 
celeration of 0.4 km/s, with a value of -25.7 km/s in 1996 and - 
25.3 km/s in 1997. Considering the uncertainties given in Table 
12 the velocity can be considered to be the same. However, a 
slight deceleration in the solar rest frame is quite possible be- 
cause of the increase in radiation pressure from the Sun be- 
tween 1996 and 1997. The LOS (line-of-sight) velocity in the 
upwind direction is the result of two antagonistic effects. First, 
selection effects linked to ionization processes increase the 
mean velocity by ionizing the slower hydrogen atoms. Second, 
radiation pressure tends to push the atoms away from the sun 
and slow them down. During most of the solar cycle, radiation 
pressure is larger than the solar gravitational attraction (Pryor 
et al. 1998). As the cycle changes from minimum to maximum, 
the interplanetary hydrogen atoms feel an increasing radiation 
pressure that becomes more efficient to slow them down in the 
upwind side of the inner heliosphere. 

The downwind velocity is mostly unchanged with a value 
equal to 21.6±1.3 km/s, which is not surprising because the hy- 
drogen distribution in the downwind side of the heliosphere is 
less affected by variations in radiation pressure. In this direc- 
tion, the hydrogen is strongly depleted by ionization effects. A 
partial filling happens beyond a few AU from the sun because 
of the relatively hot temperature of the interplanetary gas. It is 
therefore not surprising that radiation pressure plays a lesser 
role in this direction as compared to the upwind direction. 

Figure|2]shows the velocity in the solar rest frame as a func- 
tion of the angle with the upwind direction. The upwind direc- 
tion is the one defined by the isovelocity contours (Quemerais 
et al., 1999) with a longitude of 252.3°and a latitude of 8.7°. 
Changes between the two orbits are within the statistical un- 
certainties of the velocity values. 
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Table 1. Equivalent Width of the H Cells 



Year (starts in June) 
Equivalent Width (km/s) 
Equivalent Width (mA) 



1996 
5.5 ±0.1 
22.3 ± 0.4 



1997 
4.9 ± 0.1 
19.9 ± 0.4 



1999 
3.68 + 0.18 
14.9 ± 0.7 



2000 
3.61 +0.18 
14.6 ± 0.7 



2001 
3.50 + 0.12 
14.2 ± 0.5 



2002 
3.26 + 0.13 
13.2 + 0.5 



MAPS OF 1996 AND 1997 




ANGLE WITH UPWIND 




ANGLE WITH UPWIND 



Fig. 2. Variation in the LOS velocity as a function of the an- 
gle of the LOS with the upwind direction. Here, the upwind 
direction is taken to be 252.3°, 8.7°. Both curves are very sim- 
ilar showing that the 1996 and 1997 velocity profiles are very 
close and coiTespond to solar minimum conditions. The solid 
line joins the values found in 1996 and the dotted line joins the 
1997 values. 



Fig. 3. Difference in LOS velocity between the 1996 orbit and 
the five other orbits. The velocities are shown as a function of 
the angle with the upwind direction. The 1996 velocity profile 
is shown in Fig. |2l The 1997 profile (solid line, top) is very 
close to the 1996 profile. The 1999 one (solid line bottom) is an 
intermediary profile, while the 2000 (dotted line), 2001 (dashed 
line), and 2002 (dash-dot line) profiles are very similar. 



3.2. Variation from solar minimum to solar maximum 

The SWAN instrument did not operate between the end of June 

1998 and October 1998, because of the accidental loss of con- 
tact with the SOHO spacecraft. Operations resumed at the end 
of 1998 but soon were stopped after the failure of the last gy- 
roscope. Normal operations resumed in spring 1999. 

In this section, we present maps obtained from measure- 
ments made beween June 1999 and June 2003. This period 
covers the solar activity maximum of 2001 and the beginning 
of cycle 23. There are no measurements in the southern eclip- 
tic hemisphere because the -Z unit H cell stopped absorbing in 

1999 for a reason that is still unclear, while the +Z unit hydro- 
gen cell is still functional in 2006. Table^shows the variation 
in the equivalent width of the cell as a function of orbit year. 
Between 1997 and 1999, the +Z H cell shows a strong decrease 
in absorption and after that a steady but slower decrease. The 
equivalent width in 2002 is only 60% of the post launch value. 

Figure|3]shows the difference between the 1996 orbit veloc- 
ity profile and the five other orbits. The differences are shown 
as a function of the LOS angle from the upwind direction. First, 
we note that the 2000 to 2002 orbits yield very similar values, 
within 1 km/s almost everywhere. Since the 1996 and 1997 or- 
bits are also very similar, we see two main regimes for the ve- 
locity: one for solar minimum conditions with an upwind LOS 
velocity around -25.5 km/s and one for solar maximum con- 
ditions with an upwind LOS velocity around -21.5 km/s. The 
1999 and the missing 1998 ones are intermediary states. 

Table 13 gives the numerical values of the upwind velocity 
for each orbit. The values are slightly different from those in 



Fig. |3] because the averaging was done in a different way. In 
this table, velocities are averaged within 20° of the upwind di- 
rection to get better statistics. The apparent variation in upwind 
velocity can be explained partially by effects of the radiation 
pressure; however, detailed calculations will be necessary to 
see if the models match the measured variations. The general 
deceleration observed in the solar rest frame corresponds to in- 
creasing values of radiation pressure as expected from solar 
Lyman a flux measurements (Rottman et al. 2006). The change 
between 1997 and 1999 is very abrupt (3 km/s) and seems too 
important for the change of radiation pressure (h-10%) derived 
from solar flux measurements. Changes of ionization processes 
could also be partly responsible for this change of LOS ve- 
locity. Indeed, ionization processes favour fast atoms because 
slower atoms have a higher probability of being ionized. This 
heads to the selection of fast atoms and an increase in the bulk 
velocity in the inner heliosphere. However even the fast atoms 
are ionized with increasing efficiency of the ionization pro- 
cesses. This means that the atoms contributing to the intensity 
are farther away from the sun. How this wifl affect the line-of- 
sight velocity is not clear. Only time-dependent calculation of 
the hydrogen distribution will allow us to discriminate between 
the various effects involved here. 

Figure |4] presents the LOS velocity difference between 
2001 and 1996 for all directions in the northern ecliptic hemi- 
sphere. The difference in the upwind direction is roughly 4 
km/s, as seen before. The difference in the downwind direc- 
tion is small, which suggests that changes in radiation pressure 
over the cycle is not as important with the velocity distribution 
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Table 2. Upwind LOS velocity for various orbits 



Year (starts in June) 


1996 


1997 


1999 


2000 


2001 


2002 


Upwind V(year) (km/s) 


-25.7±0.2 


-25.3±0.2 


-22.5±0.5 


-2I.5±1.2 


-21.5±0.3 


-21.4±0.5 


Upwind V(year) - V(1996) (km/s) 


0. 


0.4 


3.2 


4.2 


4.2 


4.3 





SWAN/SOHO 









ECLIPTIC COORDINATES 




Fig. 4. Difference between the velocity maps of 2001 and 
1996. The maximum difference is found in the upwind direc- 
tion with a value around 4 km/s. Downwind, the difference is 
close to zero. Below 10 degrees of latitude, the results are very 
noisy. Contours of the 1996 velocity map are shifted towards 
higher latitudes than the contours of the 2001 velocity map. 
This is due to the existence of anisotropies in the hydrogen 
distribution of 1996 which create a deviation from the non- 
isotropic 2001 distribution. 



of hydrogen atoms in the downwind direction. We clearly see 
a shift between the iso-contours of velocity for 2001 and 1996. 
This is explained by the changes in the velocity contours in- 
duced by anisotropies of the solar ionization fluxes (Lallement 
et al. 2005; Koutroumpa et al. 2005). In 2001, when solar ion- 
izing fluxes are almost isotropic, the constant velocity contours 
are well fitted by cones centered on the upwind direction. In 
1996, however, due to diff'erent ionization fluxes at different 
heliographic latitudes the resulting iso-velocity contours are 
elongated towards high latitudes. This is demonstrated in Fig. 
|4]where the difl'erence of velocity maps shows a maximum that 
is shifted towards higher ecliptic latitudes. 

In this section, we have shown how the LOS velocity of the 
interplanetary Lyman a line profile changes during the solar 
cycle. The amplitude is large, more than 4 km/s in the upwind 
direction. It is also abrupt because we find a change of 3 km/s 
between 1997 and 1999. The main cause for this variation in 
the velocity is the change of radiation pressure during the so- 
lar cycle (Pryor et al. 1998). Other effects linked to ionization 
processes and their solar cycle variations may also be involved. 



4. Line-of-sight Temperatures 

This section presents the LOS kinetic temperature maps de- 
duced from the line-profile reconstruction technique presented 
in Quemerais et al. (1999). What we call the line-of-sight 



(LOS) kinetic temperature, also called apparent temperature in 
Quemerais and Izmodenov (2002), is actually the line width 
converted to temperature units using the relation between ther- 
mal velocity Vth and line width. 

Using the same notation as before, we have the foUowing 
expression 



2k k X, 



^off 



(8) 



where /(v) gives the intensity profile as a function of the LOS 
projected velocity v. Line width values are more difficult to de- 
rive from the SWAN H cell data than line shifts. Indeed, if there 
is some stellar light, noted Istar here, added to the Lyman a data 
for a given LOS, the absorption becomes 



A(measured) — l—R 



toff + ^star 



A(real) 



1 



1 + Istar/ loff 



.(9) 



If the value of Ismr is not equal to zero, the value of the 
mean line shift will not be affected much, but it will change the 
integral of the absorption profile and thus change the derived 
temperature significantly. To solve this problem, we used the 
measurements of the so-called BaF2 pixel of the +Z sensor. 
As mentioned by Bertaux et al. (1995), each detector has two 
active pixels on the sides of the 5x5 array. One of the side 
pixels of the +Z sensor was covered by a window made with 
BaF2. This window is opaque for Lyman a photons. Using the 
data recorded by this pixel over many months, we compiled a 
full-sky map that excludes the Lyman a background. This map 
can then be used to determine which areas of the sky are free of 
stellar contamination. This way, we created a mask that allows 
us to keep only the data not contaminated by starlight. 

Figure |5l shows the LOS temperature as a function of the 
upwind angle derived from the orbits of 1996 and 1997. For 
each orbit, there are 2 curves, one for angles larger than 30° 
and one for angles smaller than 30°. Because the upwind direc- 
tion is very close to the galactic plane, applying a strict limit on 
stellar counts for each LOS removes all points within 30° from 
upwind. However, some lines of sight show a low stellar con- 
tamination, of a few percent of the upwind intensity. Keeping 
those points, we were able to determine a curve for angles be- 
tween zero and 30°. This second curve shows larger uncertain- 
ties. The most striking feature of these two curves is that they 
are not monotonic. The temperature reaches a minimum value 
around 1 1000 K for an upwind angle of 60°. This feature does 
not appear in hot model calculations whether they include full 
radiative transfer effects on the line profile or simply compute 
the first-order scattering of photons (Quemerais 2000). For a 
hot model, the temperature dependence of line profiles from 
upwind to downwind is always monotonic. We also find that 
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1996 and 1997 



Solar Rest Fra 




1 Upwind (degre 




Fig. 5. LOS kinetic temperature as a function of the upwind 
angle. The values were obtained for the two orbits in 1996 and 
1997. For each orbit, there are 2 curves, one for angles smaller 
than 30° and one for angles larger than 30°. For angles smaller 
than 30°, the limit on stellar light counts has been relaxed to 
ensure that some points with a small contamination are kept. 
Contrary to hot model predictions, the temperature curve is not 
monotonic and shows a minimum value of 1 1000 K for angles 
close to 60°. The 1996 values are given by the diamonds and 
the 1997 values are given by the triangles. 



Fig. 6. Model of the line profile generated for two distinct pop- 
ulations of hydrogen atoms. One has the parameters of the in- 
terstellar gas, i.e. velocity of 30 km/s (26 km/s accelerated by 
selection effects, see text) and a temperature of 6000 K, the 
other has parameters reflecting a deceleration and heating due 
to the heliospheric interface, i.e. velocity of 20 km/s and a tem- 
perature of 14000 K. The resulting line profile in the upwind 
direction has a mean shift of -25 km/s and a temperature of a 
bit less than 14000 K. The line profiles are shown in the solar 
rest frame. 



the upwind LOS temperature is around 14000 K and the down- 
wind value is close to 18000 K. This departure from a mono- 
tonic variation with the upwind angle was also noted by Costa 
et al. (1999) using a different method. 

Quemerais (2000) computed LOS temperatures as seen 
from Earth's orbit for various models of the hydrogen distri- 
bution. What was found is that the LOS temperature of the IP 
line within 40° to 50° in the upwind direction is almost con- 
stant. The line width starts to increase for angles larger than 
50°. This increase in the line width is due to changes of the line 
shape due to acceleration and deceleration by radiation pres- 
sure and also to selection effects of the fast atoms through ion- 
ization processes. 

In that respect, we can argue that the decrease in LOS tem- 
perature between to 50° is a clear signature of the existence of 
two hydrogen populations contributing to the total line profile. 
Following Izmodenov et al. (2001), we can divide the hydrogen 
atoms in the heliosphere into four distinct populations. First, 
we consider the interstellar component that has gone through 
the heliospheric interface without any interaction with the pro- 
tons. This population at large distance from the sun has the 
distribution parameters of the interstellar gas, i.e. a bulk veloc- 
ity close to 26 km/s and a temperature close to 6000K. The 
second population is the one created by charge exchange with 
protons of the compressed interstellar plasma. Models predict 
a strong deceleration and heating of this population. These two 
populations are the main contributors to the IP line profile at 
1 AU, while the other two populations created by charge ex- 
change with the solar wind can be neglected here (Quemerais 
and Izmodenov 2002). 

A simple model of the two populations is shown in Fig.|6] 
The profiles are represented by Gaussian functions. The plots 
are made in the solar rest frame. One population has the pa- 



rameters of the interstellar gas, i.e. a velocity of 30 km/s and 
a temperature of 6000 K. We used a velocity of 30 km/s for 
the primary component to account for selection of faster atoms 
by ionization processes (See Quemerais and Izmodenov 2002, 
Table 4). The other component is decelerated in the solar rest 
frame and heated. Its parameters are given by a velocity of 20 
km/s and a temperature of 14000 K. We have computed the line 
shift and line width of the sum of these two line profiles pro- 
jected on an LOS with an angle from upwind between zero and 
50°. The results are shown in Table|3l We find values similar to 
the observations. 

We do not claim that this simple model fits the data. It is just 
an example to illustrate why the LOS temperature decreases 
between and 50°, i.e. because the Doppler shift between the 
two components of the line decreases as the cosine of the angle 
with the upwind direction. After 50° from upwind, dynamic ef- 
fects on the hydrogen distribution make the line width increase 
again. Actual modeling is required here to correctly interpret 
these data. 

Figure shows the LOS temperatures found for the four 
orbits from 1999 to 2002. The data from these orbits are nois- 
ier due to the degradation of the sensitivity of the sensor units 
(Quemerais and Bertaux 2002). It was not possible to recover 
profiles for LOS's within 20° of the upwind direction or within 
30° of the downwind direction. However, two results can be 
seen from Fig.0 First, the temperature minimum around 60° 
from upwind has either shifted to lower angle values or even 
disappeared. This suggests that the shift between the two com- 
ponents is smaller than in 1996 or that one of the components 
has become less important relative to the other giving more 
weight to the parameters of the other component in the line 
profile. Second, the curves show lower temperatures than what 
was found in 1996 and 1997. The overall decrease is around 
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Table 3. LOS velocity and temperature for a 2-population model 



LOS Angle from Upwind 


LOS Velocity 


LOS Temperature 


0° 


-23.9 km/s 


13733 K 


10° 


-23.6 km/s 


13646 K 


20° 


-22.5 km/s 


13394 K 


30° 


-20.7 km/s 


13008 K 


40° 


-18.3 km/s 


12535 K 


50° 


-15.4 km/s 


12032 K 


60° 


-12.0 km/s 


11559 K 



LOS Temperature 
2.2x10' I ' ' ' ' \ ' ' ' ' \ ' ' 




50 100 150 200 

ANGLE LROM UPWIND in degrees 



rower line profiles. These results will be compared with model 
computations in a future work. 

5. Comparison with HST profiles 

Direct measurements of the Lyman a line profile require a 
very good spectral resolution which is not easily obtained with 
a space instrument. Fortunately, the STIS instrument on the 
Hubble Space Telescope has provided a few measurements in 
June 2000 and March 200 L This section presents the available 
upwind spectrum and a comparison with the results obtained 
from the SWAN H cell data. 



Fig. 7. LOS Temperature as a function of the upwind angle. 
The values were obtained for the four orbits from 1999 to 2002. 
The data from 1999 are shown by crosses, those of 2000 by 
stars, those of 2001 by squares, and the data of 2002 are shown 
by triangles. There are no values for angles smaller than 20° 
or larger than 150° because of contamination by stellar light. 
The two curves from 1996 (dashed line) and 1997 (dotted line) 
have been added for comparison. First, the temperature mini- 
mum seen previously around 60° seems to have shifted towards 
lower values. Second, the profiles appear cooler by roughly 
1000 K everywhere. 



1000 K. This feature also suggests that one of the components 
has become relatively less important, thus yielding an appar- 
ently cooler profile. It is unlikely that the shift between the 
components has decreased, because this wouldn't change the 
LOS temperature for crosswind lines of sight. 

The LOS temperature profiles presented in this section 
show that the IP line profile is most likely made up of two 
components due to distinct hydrogen populations. The first one 
is composed of unperturbed interstellar hydrogen atoms get- 
ting close to the sun. The second one is composed of hydrogen 
atoms created after charge exchange with interstellar protons 
compressed and heated in the heliospheric interface. The bulk 
velocity difference explains why the IP line profile line width 
decreases when the upwind angle of the LOS goes from zero 
to 50°. We also find that the profiles appear cooler during solar 
maximum than during solar minimum. This could indicate that 
the increase of the ionization processes but also radiation pres- 
sure with the solar cycle is more effective on one population 
than the other The slower population will be more ionized than 
the fast one for instance. Consequently, this will result in nar- 



5.1. The data 

The data presented in this section were obtained with the STIS 
instrument on-board the Hubble Space Telescope. One mea- 
surement was obtained in March 2001 (Upwind LOS). 

The main problem for the observation of the interplane- 
tary line profile from Earth's orbit is caused by the existence 
of the strong emission of the geocorona. At the altitude of the 
Hubble Space Telescope, the geocoronal emission is 5 to 15 
times brighter than the interplanetary line depending on the di- 
rection of the LOS. The best time to observe the interplanetary 
line is when the Doppler-shift between the two lines is at its 
maximum. For the upwind direction, this is in March when the 
Earth velocity vector is toward the upwind direction. In that 
case, the relative motion between the H atoms and the observer 
is close to 50 km/s. 

Crosswind observations are limited to a relative velocity of 
30 km/s because their LOS is perpendicular to the velocity vec- 
tor of the interstellar H atoms. This maximum value is obtained 
when the Earth velocity vector is perpendicular to the interstel- 
lar wind direction, i.e. either when the Earth is upwind from 
the Sun (early June) or downwind (early December each year). 

The upwind measurement was obtained on March 21, 2001 
when the Doppler shift between the interplanetary line and the 
geocorona is close to its maximum (Position 2). In that case, the 
two lines are well separated. Figure|Slshows both the geocorona 
and the interplanetary line. A symmetric line shape fitting the 
geocorona has been removed from the data. The residual shows 
the upwind interplanetary line. 

We estimated the LOS values of the velocity and the tem- 
perature of the line. The results for apparent velocity corrected 
for the Earth's motion is 

Vlos = -20.4 + 0.2 km/s , 
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Fig. 8. Upwind Interplanetary Spectrum measured by STIS on 
March, 2 1 200 1 . The geocorona (dash line) is more intense than 
the IP line (factor of 6) but the two lines are well separated be- 
cause of the Doppler shift between the H atoms and the Earth. 
The geocoronal emission has been fitted and removed, leaving 
a larger uncertainty longward of the interplanetary line. 

the LOS temperature is 

Tlos = 16570 + 500 K. 

Those values are larger than the actual values because of the 
convolution of the actual line profile with the line spread func- 
tion of the instrument. The geocoronal emission from Fig. |S] 
gives a good estimate of the line spread function (LSF). Its LOS 
temperature is equal to 5300 K whereas the actual temperature 
of the geocorona is around 1000 K. 

We have used the LSF deduced from the geocorona to de- 
convolve the upwind Une profile. The result is shown in Fig.|9l 
First, the data were fitted to a Voigt function. Then, this func- 
tion was deconvolved yielding the spectrum. By assuming that 
the geocorona gives the LSF, we have slightly overestimated 
the actual width of the LSF, as reflected in the larger uncertainty 
in the temperature estimate. Future observations of the martian 
neutral H atom emission at Lyman a will better estimate the 
LSF because the martian emission profile has a thermal width 
equivalent to a few hundred K (« 200 K). 

The resulting parameters for the upwind LOS are after de- 
convolution, 

y^os = -20.3 + 0.2 km/s . 

The LOS temperature is 

r^os = 10970 + 1000 K . 

The results found from the HST upwind line profile are com- 
patible with the SWAN H cell measurements. First, the mean 
line shift measured in 2001 was 20.3 km/s. The SWAN result 
is 21.4 km/s. Taking a possible bias into account due to the re- 
moval of the geocoronal line from the HST spectrum, we get a 
correct agreement, thus confirming the change in the line shift 
of the IP line from solar minimum to solar maximum. Note also 
that a previous HST measurement made by GHRS (Clarke et 
al 1995, 1998) agreed with the SWAN value of -25.7 km/s for 
the solar minimum IP mean line shift. 



Fig. 9. Upwind interplanetary spectrum in the solar rest frame. 
The diamonds show the data, the thin solid line shows a Voigt 
function fit to the data, the thick solid line shows the decon- 
volved spectrum obtained from the Voigt function fit assuming 
that the LSF is given by the geocoronal profile. 

The temperature found from the profile is close to 1 1000 K. 
This value may be slightly underestimated because the LSF of 
the STIS instrument is not as wide as the Earth's coronal line. 
Comparing with the values shown in Fig. we find correct 
agreement. This also confirms that the upwind IP line profile 
LOS temperature has decreased from 14000K to 1 lOOOK from 
solar minimum to solar maximum It also shows that the tem- 
perature inflexion seen at 60° from upwind at solar minimum 
more or less disappeared at solar maximum. 

6. Conclusion and discussion 

This analysis of the SWAN H cell data has allowed us to recon- 
struct interplanetary Lyman a line profiles between 1996 and 
2003. This period covers the solar activity minimum of 1996 
and the maximum of 2001. 

We have found that the mean line shift changes from a LOS 
velocity of 25.7 km/s to 21.4 km/s in the solar rest frame. This 
deceleration is mainly due to changes in radiation pressure with 
increasing activity, although changes of the ionizing fluxes are 
also involved. Detailed modelling will be necessary to repro- 
duce this large variation in the mean line shift. 

A comparison with a spectrum recorded by STIS on HST 
yields a good agreement. The STIS line mean shift corresponds 
to an LOS velocity of 20.3 km/s which, given the uncertainties 
and possible biases involved in both analyses, is quite accept- 
able. We should also point out that the mean line shift change 
seen by SWAN is very rapid. We find a variation by 3 km/s 
between 1997 and 1999. This seems hard to explain only by 
changes in radiation pressure. 

The SWAN H cell data were also used to determine line 
widths (LOS temperature). It is found that the line width vari- 
ation with upwind angle is not monotonic as is usually found 
from hot model computations. This can be explained as a proof 
that the IP line profile is made of two distinct components scat- 
tered by populations with different bulk velocities and tem- 
peratures. These different components have been theoretically 
predicted by models of the hydrogen interaction of the heHo- 



QUEMERAIS ET AL.: IP LINE PROFILES 



9 



spheric interi'ace. Here we have an observable effect on the Une 
width which is created by these two populations. Actual model 
computations of hydrogen distribution and backscattered pro- 
files will be made to test this explanation. 

Finally, we found that the LOS temperature profiles are 
cooler during solar maximum than solar minimum. The line 
width change corresponds to a decrease in the LOS tempera- 
ture by 1000 K. A tentative explanation is that the slow hy- 
drogen population component is more effectively ionized than 
the fast one during solar maximum. This results in a smaller 
contribution to the total Une shape and hence a narrower Une 
profile. 

The results obtained in this analysis are sunmiarized by this 

Ust 

- Between 1996 to 2001, the mean line shift of the interplan- 
etary Lyman a line changes from a LOS velocity of 25.7 
km/s to 21.4 km/s in the solar rest frame. 

- The Une shift changes abruptly between 1997 and 1999 by 
3 km/s. 

- During solar minimum, the IP Lyman a line width shows a 
variation with the angle from upwind, which is not mono- 
tonic but has a minimum around 60°. This suggests that the 
IP line is composed of two components with different mean 
line shifts. 

- During solar maximum, the LOS kinetic temperatures de- 
crease slightly and the minimum is less pronounced, sug- 
gesting that the ratio between the different components is 
changed from solar minimum conditions. 

- Spectra obtained by HST in 1995 and 2001 give LOS ve- 
locities and temperatures are compatible with the SWAN H 
cell measurements. 

These empirical results wiU be confronted to model calcu- 
lations in a future work. 
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